Symmetries and infinitesimal symmetries 
of singular differential equations 



Xavier Gracia a and Josep M. Pons 6 

a Departament de Matematica Aplicada IV 
Universitat Politecnica de Catalunya 
Campus Nord UPC, edifici C3 
C. Jordi Girona, 1 
08034 Barcelona 
Catalonia, Spain 

b Departament d'Estructura i Constituents de la Materia 
Universitat de Barcelona 
Av. Diagonal, 647 
08028 Barcelona 
Catalonia, Spain 

e-mails: xgraciaOmat. upc.es, pons@ecm.ub.es 
15 January 2002 

Abstract 

The aim of this paper is to study symmetries of linearly singular differential equations, 
namely, equations that can not be written in normal form because the derivatives are 
multiplied by a singular linear operator. 

The concept of geometric symmetry of a linearly singular differential equation is intro- 
duced as a transformation that preserves the geometric data that define the problem. 
It is proved that such symmetries are essentially equivalent to dynamic symmetries, 
that is, transformations mapping solutions into solutions. Similar results are given for 
infinitesimal symmetries. 

To study the invariance of several objects under the flows of vector fields, a careful 
study of infinitesimal variations is performed, with a special emphasis on infinitesimal 
vector bundle automorphisms. 
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1 Introduction 

Let M be a manifold, and V the set of paths in M. Consider a subset S C V, which 
may be understood as the set of solutions of a certain problem on the set of paths. This 
problem — usually we think of a differential equation — may be stated in terms of several 
objects, the data of the problem. The solutions of the problem are the paths in V that 
satisfy some condition with respect to the data that identify the problem. Examples of 
data can be a vector field, a submanifold, a connection, a potential function, a lagrangian, 
etc, and for each case a specification of the associated problem must be given: let it be the 
search for the integral curves of a vector field, the critical paths of an action functional, 
etc. 

Since the problem is identified by the data and solved by the solutions, there appear 
naturally two different concepts of what a symmetry of the problem is, differing as to 
whether the emphasis is put on the side of the data or on the side of the solutions. In 
a certain sense, a transformation that preserves the data is a symmetry of the problem; 
and so it is, but in another sense, a transformation that maps solutions into solutions. In 
order to avoid confusion, we can call geometric the former symmetries, and dynamic the 
latter. In fact, the word dynamic has the reminiscence of the equations of motion that set 
the dynamics of a physical problem. Also, in most cases, the data of the problem will have 
geometric significance, and hence the name suggested. Usually a geometric symmetry of 
a problem will also be a dynamic symmetry; thus the search for geometric symmetries 
will be a relevant part of a wider subject: the search for dynamic symmetries. Noether 
transformations of an action functional are an example of this. 

As for the type of transformations of the paths, we will confine ourselves to point 
transformations (of the dependent variable), which arise as 7 1— > 99* [7] := ^07, for a certain 
diffeomorphism <£>: M — > M. (Among non-point transformations we have for instance 
reparametrisations of the independent variable, and the generalised transformations, where 
the transformation involves also the derivatives of the path — see | 01v93(| .) 



In this paper we will deal with problems resulting in first-order autonomous differen- 
tial equations on M. Among the possible data that identify the problem, there is the 
differential equation itself, considered as an implicit relation involving a path x(t) and its 
derivative with respect to the evolution parametre: 

F(x,x) = 0. 

Equations of this form are often called differential-algebraic, or implicit differential equa- 
tions. 

Of course, if one can isolate the derivative, 

x = f(x), 

the equation is said to be in normal form, and giving an initial condition (to,x(to)) de- 
termines uniquely the solution x(t). However, we are mainly interested in the singular 
case. 
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More precisely, we are interested in implicit differential equations of the form 

A(x)x = b(x), 

where the velocities can not be isolated because of an everywhere singular linear oper- 
ator A{x) multiplying them. Such equations may be called linearly singular differential 
equations. This general class of implicit differential equations was geometrically presented 
in [ GP 9 1 1 |GP 92(| . In these papers it is pointed out that many interesting systems of 
mathematical physics and applied mathematics are linearly singular. 

In more detail, this framework describes the equations of motion of the presymplec- 
tic dynamical systems [ GNH 78}] (including their applications to lagrangian and hamilto- 
nian mechanics [Dir64] [GN 7£] [ MT 78| ] |Ski83 | [ [5R83] ]), the first-order lagrangian formal- 
ism [GP92], the higher-order lagrangian dynamics |GPR91 | [ LR85[| (including also their 
"higher order differential equation" conditions [GPR92]), and systems with nonholonomic 
constraints. 

In addition to these applications of interest for mathematical physics, one can find ap- 
plications of implicit and linearly singular differential equations to electrical and chemical 
engineering, control theory, economics, etc (see for more details examples and references in 
[ GP 92| pMR 96f §ILR8|] [|Rhe 84|| ). As well as all the mentioned papers, there are also 
many articles and books studying geometric features |CQ88fl []HB 84fi [fvIMT 92fl RMMT 95fl 
[ MR 991 [|MT T8|1 [|Rei 9C|| [p,ei91| |[RR94jl [|Tak 76| and numerical methods [|Cam 8§ [ |HLR89|1 
[HW91] for implicit equations. 

The symmetries of an implicit differential equation can be studied using general tech- 
niques [01v93]. Important topics as lagrangian systems and presymplectic systems have 
been widely studied (see for instance | 01v93 | | LM 96| and references therein). Besides 



these cases, there are few references on symmetries of implicit differential equations: we 
could point out the articles [ MMT 92| , where symmetries and constants of motion for im- 
plicit systems F(x,x) = are studied; [pO 88 |, which contains a study of normal forms 
of linearly singular systems given by a vector bundle morphism A: TM — > TM; [MR9S], 
which deals with symmetries of linearly singular systems given by a vector bundle mor- 
phism A: TM — > T*M; and the recent study [BS 01 1 of symmetries and reduction of Dirac 
structures. 

The main purpose of this article is to study the symmetries of a linearly singular 
differential equation. For such an equation we can consider the geometric symmetries 
preserving the data (A, b) that define the equation. It will be proved that any dynamic 
symmetry of the differential equation may be locally realised as a geometric symmetry of 
the data. A similar result will be also given for infinitesimal symmetries. 

The paper is organised as follows. Section 2 presents the geometric framework of 
linearly singular differential equations, and gives several basic results. Section 3 studies the 
symmetries of such a system, and relates them to the symmetries of the associated implicit 
differential equation. In section 4 the concept of infinitesimal symmetry is presented, and 
a study similar to that of section 3 is performed. Section 5 particularises all the results to 
regular and to consistent systems. Section 6 is devoted to an example, and section 7 to 
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conclusions. Finally, there is an appendix dealing with calculus of infinitesimal variations, 
and more particularly to the invariance of maps under the action of flows of vector fields, 
and to infinitesimal symmetries of vector bundles. 

The tools used in this paper are those of differential geometry, in particular manifolds 
and submanifolds, vector fields and their flows, and vector bundles and their morphisms 



| AMR 83 ] [ Die 7C | KMS 93 1. Throughout the paper the manifolds are finite-dimensional 
and paracompact, and the maps are smooth; "differential equation" means "first-order 
autonomous ordinary differential equation" . 

2 Linearly singular differential equations 



In this section we recall some of the main results from [ GP 91 [ GP 92 ], and we give 
additional results to be used later on. 

2.1 The geometric framework 

Definition An implicit system on a manifold M is a submanifold D C TM. It defines 
an implicit differential equation, for which a path / —* M is a solution when its lift to 
the tangent bundle, £, is contained in D: 

e(J) C D. (2.1) 

When D is the image of a vector field X on M, D = X(M), one has an explicit 
differential equation, or says that the equation can be put in normal form. Then a path £ 
is a solution of D iff £ = X o £. 



Definition ([ |GP 91 ] [ GP 92[ ]) A linearly singular system is a quintuple (M, F, n, A, b) 



given by a manifold M, a vector bundle ir: F — > M, a vector bundle morphism A: TM — > F, 
and a section b: M — > F. These data define a linearly singular differential equation, for 
which a solution is a path £: J — > M such that 

Ao^ = 6o£. (2.2) 

The associated implicit system is the subset 

D = A- 1 (f,(M)) C TM. (2.3) 

We will use the notation (A: TM — ► i 7 , 6) to refer to a linearly singular system. The 
following diagram shows all these data: 

TM — F 
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Proposition 1 The differential equations defined by a linearly singular system 
(A: TM — ► F, b) and its associated implicit system D have the same solutions. 

Proof. Equation (2.1) means that, for each t, £(t) £ D. This is equivalent to A-£(t) € 
b(M), and being A fibre-preserving this is equivalent to A-£(t) = &(£(£)). ■ 

Note that in general ( |2.2| ) may not have solutions passing through every point in M, 
and if there is a solution passing through a point x at a given time it may not be unique. 
We call motion set the set 5cMof points by which a solution passes. 

It is useful to try to describe the solutions of the equation of motion ( p^ ) as integral 
curves of vector fields. More precisely, if M 1 C M is a sub manifold and X is a vector field 
on M tangent to M', then the integral curves of X contained in M' are solutions of the 
equation of motion (2.2) if and only if X satisfies 

AoX ~ b, (2.4) 

M' 

where the notation means equality at the points of M'. Let us remark that this is an 
equation both for X and M', since in general there will not be a vector field satisfying 
this equation all over M. 



2.2 The constraint algorithm 

Consider a linearly singular system (A:TM — > F,b). To solve the corresponding differ- 
ential equation a consistency algorithm may be performed. This algorithm is indeed a 
generalisation of the presymplectic constraint algorithm |GNH 78 ], which is a geometri- 
sation of Dirac's theory for singular lagrangians [ Dir 64 j. Let us describe this algorithm 
briefly. 



Definition The primary constraint subset is the set M\ C M of points x where the 
linear equation A x -u x = b(x) is consistent: 

Mi = {x G M | b(x) € ImA x }, (2.5) 

The functions of M vanishing on Mi constitute the ideal of primary constraints. 

The reason for the terminology is clear: in view of the differential equation of motion 
( |2.2| ), if £ is a (smooth) solution of it, then necessarily £ lives in M\. As for the constraints, 
in this paper we are not especially interested in explicit procedures to compute them; see 
[ GP 92| for more details. 

To proceed further it is convenient to require some regularity conditions on A and Mi : 

Definition We will refer to the regularity assumption as the following conditions to be 
satisfied by a linearly singular system (A: TM — > F, b): 

1. The morphism A has constant rank — thus Ker A C TM and ImA C F are vector 
subbundles, D C TM is a closed submanifold and Mi is a closed subset. 

2. The primary constraint subset Mi C M is a non-empty submanifold. 
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Let us assume that our linearly singular system satisfies the regularity assumption, 
and let £ be a solution of the corresponding differential equation. It is clear that £ is a 
also a solution of the linearly singular differential equation defined by restricting all the 
problem to M\, namely, the subsystem (A\:TM\ — * Fi,bi), where F\ = F\ Mi , and A\ 
and b\ are the corresponding restrictions to M\. 

Note that the problem is not yet solved: for a point x € Mi, b(x) does not necessarily 
belong to the image of the restriction of A x to the subspace T^Mi C T X M. Repeating 
the consistency analysis for the subsystem yields a subset M2 := (Mi)i C M\. 

Let us assume that the regularity assumption holds for the successive subsystems. The 
repetition of the consistency analysis on the subsystems yields an algorithm that reaches 
— in a finite number of steps since M is finite-dimensional — a final constraint submanifold 
Mf := P| Mj. The solutions of the original problem are the solutions of the equation of 

motion of the linearly singular system (Af : TMf — > Ff , bf) defined by restriction to Mf . By 
construction, bf has its image in ImAf — otherwise the algorithm would not be finished. 
Therefore the equation 

A f oX f = b f (2.6) 

for a vector field Xf in Mf has solutions. Since Mf is closed these solutions can be extended 
throughout M to yield solutions X of the equation of motion (2.4) along M' = Mf which 
are tangent to this submanifold. Given a particular solution Xf of ( p.6[ ), the set of solutions 
is Xf + Ker Af. Therefore there is a unique solution (on Mf) iff Af is injective. 



TMf 



Im A { 



Ff 



X 




M f 

So note that the final dynamics is simply that of a linearly singular system where the 
morphism A is surjective. 

Note finally that if the Mj fail to be submanifolds, then in order to apply the constraint 
algorithm some points of the base space M may have to be removed; in this case only a 
subset Mf C S of the motion set will be obtained, and the motion set S may not be a 
submanifold — see some examples in [GP92]. 



2.3 Morphisms of linearly singular systems 

Definition A morphism of linearly singular systems between (A: TM — > F, b) and 
(A': TM' — ► F' , b') is a morphism (</?, <£) between the vector bundles F — > M and F' — > M' 
(so it satisfies o 7T = ir' o <E») such that 

$06 = 6' 099, (2.7.a) 
$»i = A'oTip. (2.7.b) 
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Let us show all this in a diagram: 



TM 




M 

With this definition the linearly singular systems constitute a category. Its isomor- 
phisms correspond to the case when (<p, <I>) is an isomorphism of vector bundles. In this 
case, in general we can define 



$^[A] := $oio (TV)" 1 , $*[&] := $o6of 



(2.1 



where <3?* denotes the push- forward through the isomorphism (<p, $); then the condition 
to be an isomorphism of linearly singular systems is 3?*L4] = A', <]?*[£>] = b' . 

Note also the following trivial remark: if <&: F — > F is a base-preserving automorphism, 
it defines an isomorphism between (A: TM —* F, b) and (<&oA: TM — > F, $06). This reflects 
the fact that the equations Ao£ = 6o£ and (^oA)o^ = ($o&)o£ are completely equivalent. 

It is easily proved that a morphism applies solutions of the corresponding differential 
equation into solutions. Other constructions with linearly singular systems can be carried 
out: subsystems, quotients, products, . . . These constructions induce natural morphisms. 



See |GP 92|| for more details. 



2.4 Primary dynamical vector fields 

Let us have a closer look to the first stage of the constraint algorithm. 
Proposition 2 Consider a linearly singular system {A:TM — > F,b). Then: 

1. M\ = tm{D), where tm- TM — > M is the natural projection. 

2. If the regularity assumption is satisfied, D C TM\ M — ► M\ is an affine subbundle 
modelled on (Ker^4)| Mi . 

Proof. For the first assertion, x £ Mi iff there exists v x € T X M such that A-v x = b(x), 
which is equivalent to saying that v x £ D x — we write as usual D x = D n T X M. 

Now consider the restriction of A to the submanifold Mi, Ai: TM\ M — > F\ M . Since 
the section bi = b\ M is in the image of Ai and Ai is a vector bundle morphism with 



constant rank, we have that D 
on Keryli = (KerA)\ M . 



A 1 (61 (Mi)) is an affine subbundle of TM\ M modelled 
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Definition A section of D — > M\ is called a primary dynamical vector field. 

Of course we can suppose that such a section is extended to a vector field X on M. 
Then, saying that X is a primary vector field means that 

AoX ~ b. (2.9) 

Mi 

Such vector fields constitute a first approach to the final dynamics (the tangency of X to 
Mi is not guaranteed). 

If X a is a primary dynamical vector field, another vector field X is a primary field 
if and only if it differs from X Q on a section of Ker^4 on M\. Thus, if (r„)i<„< m is a 
local frame for Ker A near M±, then there locally exist functions g^, uniquely determined 
on Mi, such that, locally, 

X t, ^o + E^; (2.10) 



see [GP92] for more details on how an explicit computation of the final dynamics can be 
obtained in this way. 



3 Symmetries of linearly singular systems 

We have pointed out in the introduction that one may define several concepts of symmetry 
of a differential equation according to the data that define it. Our purpose now is to study 
the natural symmetries of linearly singular systems. 

Definition A symmetry of an implicit system D C TM is a diffeomorphism p: M — ► M 
leaving D invariant, that is, 

(T<p)(D) C D. (3.1) 



Proposition 3 A symmetry ip of D maps solutions of the corresponding differential equa- 
tion into solutions. 

Proof. It is immediate: if £ is a solution (that is, £(t) £ D) then ip o £ is also, since 
(<poO-(t) = (T<p)-£(t)eD. m 

Though this is the natural geometric definition of a symmetry of D, this is not a 
necessary condition for cp to define a symmetry of the solutions of the differential equation, 
because an implicit differential equation may not have solutions passing through every 
point in D. But with a convenient refinement of it, this condition essentially characterises 
the symmetries of the solutions of the differential equation. More precisely, following the 
terminology of [01y_93], let us call D locally solvable if for each v S D there is a solution 
£ of the implicit differential equation such that £(0) = v (see also | MMT 92fl ). 

Proposition 4 Suppose that D is locally solvable. Then a diffeomorphism p: M — > M is a 
(geometric) symmetry of D iff it is a (dynamic) symmetry of the corresponding differential 
equation. 
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Proof. The direct implication is the preceding proposition. For the converse, let v £ D. 
For a certain solution £, £(0) = v, and since fo£ is also a solution, (Tip)-v = (<p°£)'(0) € D. 

m 

When D is not locally solvable, if one can perform a "constraint algorithm" to pass to a 
locally solvable D' C D, then the dynamic symmetries of D (or D') are in correspondence 
with the geometric symmetries of D' . 

Definition A symmetry of a linearly singular system (A: TM — ► F, b) is an isomorphism 
with itself, that is, a vector bundle automorphism (if, <J>) of 7r: F — » M such that 

6 = $46], A = $*[A]. 



We have already said that such a transformation maps solutions into solutions. A more 
precise result is: 

Proposition 5 Let f be the base map of a symmetry (f, <&) of a linearly singular system. 
Then f is a symmetry of the associated implicit system D. 

Proof. We have to show that D is c/2-invariant. Let v x S D: A x -v x = b(x). Then, 
according to (p77|) , 

\(x)-(T x (<f)-v x ) = § x -(A x -v x ) = $ x -b(x) = b(if(x)), 

which shows that T x (f)-v x € D. ■ 

We want to prove a kind of converse of this statement. To this end, first we state an 
auxiliary result: 

Proposition 6 Let A: E — > F be a vector B-bundle morphism, and A':E' — > F' a vector 
B' -bundle morphism. Let S: E — > E' be a vector bundle isomorphism over a map ip:B — » 
B' , and such that S-KerA C Kei A' . Suppose that A and A' have the same constant rank, 
and that F and F' have the same rank. 

Then locally there exists a vector bundle isomorphism T:F^F' such that A' o S = 
To A. 

So this proposition deals with the commutativity of the "upper" square of the following 
diagram by means of a certain morphism T: 




Proof. Since S:E — > E' is an isomorphism, the hypotheses on the kernels and on the ranks 
imply that S{Kei A) = KerA', and therefore S defines an isomorphism S: Ej Ker A — > 
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E' / Ker A'. Using the canonical isomorphisms Ej Ker A = Im A, S defines an isomorphism 
T :Imi -> ImA'. 

Since any vector subbundle is a direct factor, we can put F = ImA © F\ and F' = 
Im A' © F[ . Any vector bundle map T\ : F\ — > F[ can be combined with To to obtain a 
vector bundle morphism T such that T o A = A' o S. The condition on the ranks locally 
allows to choose T\ to be an isomorphism, and therefore T. ■ 

Note however that F\ and F[ in the proof need not be isomorphic, even if F and F' 
are so. Therefore the last assertion in the proposition is necessarily local. The proof also 
shows that T is uniquely defined only on ImA 

Now we are ready for the main result: 

Theorem 1 Let (A: TM — > F, b) be a linearly singular system satisfying the regularity 
assumption, and let D = A^ 1 (b(M)) be the associated implicit system. 

Let ip: M — > M be a diffeomorphism. The following statements are equivalent: 

1. p is a symmetry of the implicit system D C TM. 

2. The restriction o/Ti^(Ker ^4) to M\ is in Ker A and, for any primary vector field X, 
the restriction of ip*[X] — X to M\ is in Ker A. 

3. cp is locally the base map of a symmetry of the linearly singular system. 
Proof. Condition (1) means that 

T x (<p)-D x = D v(x) (3.2) 
for each x € M\. Thanks to the affine structure of D — > M\ (proposition^), 

d \mi = X \m x + Ke^U' 
where X is a primary dynamical vector field on M. At x we have 

D x = X(x) + Ker A x . 

So by condition (1) we have 

T x ((p) ■ Ker A x = Ker A^ (3.3.a) 

T x (ip) ■ X(x) - X(<p(x)) £ Ker A^) (3.3.b) 

thus obtaining condition (2), and conversely. 

Now let us apply proposition ^ to the following diagram: 
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We conclude the local existence of a vector bundle isomorphism over tp\ M : M\ — ► M\ 
closing the "upper square" in the diagram. We extend &\ Ml locally to a vector bundle 
morphism <£: F — > F, which we can assure to be an isomorphism at least on an open 
neighbourhood of M\ C M. We have 

$ x oA x = A ip{x) oT x (ip), (3.4.a) 

which means that $*[j4] = A. Moreover, 

$ x • b(x) = $ x ■ A x ■ X(x) = A v{x) ■ T x (ip) ■ X(x) = A v(x) ■ X{tp{x)) = b{(p(x)), 

so we have 

$ x • b(x) = b(cp(x)), (3.4.b) 

which means that 3>*[i>] = b. 

Conversely, from = A we have the first of equations ( |3.3| ), and, from $*[6] = b, 

A^ {x) ■ (T x ((p) ■ X(x) - X(<p(x))) =$ X -A X - X(x) - A^ {x) ■ X(<p(x)) =<£> x -b x - 6„ (s) = 0, 

therefore we have the second of equations (^) . ■ 
A final remark: as before, $ is uniquely defined on Im-A^. 



4 Infinitesimal symmetries 

The infinitesimal version of an automorphism of a differential manifold is a vector field X, 
in the sense that integration of it yields a local 1-parametre group of diffeomorphisms, 
the flow Fx- We can translate the results of the preceding section to the infinitesimal 
language. The basic geometric tools are gathered in the appendix. 

Definition An infinitesimal symmetry of an implicit system D C TM is a vector field 
V on M such that the maps Fy are locally symmetries of D. 

Proposition 7 V is an infinitesimal symmetry of D iff its canonical lift to TM, V T , is 
tangent to D. 

Proof. It follows from the definition of infinitesimal symmetry and from the definition of 
the vector field V T , whose flow is constituted by the maps T(Fy). ■ 

Definition An infinitesimal symmetry of a linearly singular system (A: TM — > F, b) is 
an infinitesimal automorphism (V, W) of the vector bundle ir: F — > M such that its flow 
(Fy,F|^) is constituted by local symmetries of the linearly singular system. 

According to proposition [l4| in the appendix, being (V, W) an infinitesimal automor- 
phism of vector bundles means that V is a vector field on M, W is a vector field on F, 
and W projects to V and is a linear vector field. 

Theorem 2 An infinitesimal vector bundle automorphism (V, W) is an infinitesimal sym- 
metry of the system iff 

TboV = Wob, (4.1) 
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TAoV T = W oA. (4.2) 



Proof. The conditions (2.7) for the couple of flows (Fy,Ff^) to be a symmetry may be 
written as 



According to proposition |12| in the appendix, these equalities hold for each e iff ( |4.lD and 
(U) also do. ■ 

Now we can state the infinitesimal version of theorem [l]: 

Theorem 3 Let (A: TM — > F, b) be a linearly singular system satisfying the regularity 
assumption, and let D = A^ 1 (6(M)) be the associated implicit system. 

A vector field V on M is an infinitesimal symmetry of D iff, locally, there exists a 
vector field W on F such that (V, W) is an infinitesimal symmetry of the linearly singular 
system. 

Proof. Consider local coordinates (x l ) on M, {x l ,u l ) on TM and (x l ,v k ) on F. Then 
the section b reads (x l ) i— > (x l ,b k (x)), and the morphism A reads (x l ,u l ) i— > (x*, A^x)^). 
Let us write 

V - a i - 

dx l ' 

so that V T = a 1 -—^ + — — tV? — —. The subset D C TM is locally defined by the vanishing 

ox 1 ax- 5 ou 1 

of the constraints tp k (x,u) := A\{x)u % — b k (x). In this way, the tangency of V T to D is 
locally expressed as 

V T -^ k = Bty 1 , 

for some functions B k ; in principle, these functions depend on (x*,u 4 ), and may not be 
unique. However, the derivative of an affine function with respect to a linear vector field is 
again an affine function, so the functions B k can be assumed not to depend on u l . Writing 
more explicitly the preceding equality we obtain 

Ak^ a \.i i • j db k j k , j k ,l 

■Aa— — tU + — ra J U — ra =DiA;U —0,0. 

3 dx l dxl dx l 11 1 

Equating the constant and the linear parts, we conclude that V is an infinitesimal sym- 
metry of D iff there are functions B k (x) such that 

On the other hand, a linear vector field W on F projecting to V is expressed as 

for some other functions B k . Then, according to theorem ^, the conditions of being (V, W) 
an infinitesimal symmetry of (A: TM — » i 7 , b) read as follows: T6 o V = W o 6 means 
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dxi 3 dx l / ' 1 



and TA o V T = W o A means 



Comparing these conditions with ( [4.3|) , we conclude that being V an infinitesimal symme- 
try of D is equivalent to the existence of W making (V, W) an infinitesimal symmetry of 
(A:TM -> F,b). m 



5 Regular systems and consistent systems 
5.1 Regular systems 

Definition A linearly singular system (A: TM — > F, b) is regular if A is a vector bundle 
isomorphism. 

In this case the dynamics is uniquely determined by the associated explicit system 
given by the vector field A = A^ 1 o b. 

Note that, for a diffeomorphism <p: M — > M, now it is equivalent to say that D = A(M) 
is invariant by Tip, and that X is invariant by ip. 

If (ip, is a symmetry of the linearly singular system, then $ is uniquely determined 
from ip as 

$ = AoT^oA" 1 . (5.1) 

Then the relations $ o b = b o 93 and Ty? oI = Ioj3 are readily seen to be equivalent. So 
we have proved the following result: 

Proposition 8 Suppose that the system (A: TM —* F, b) is regular, and let p:M^M 
be a diffeomorphism. Then the following statements are equivalent: 

1. ip is a symmetry of the associated implicit system D. 

2. ip leaves the dynamical vector field X invariant. 

3. ip is the base map of a symmetry {ip, <£) of the linearly singular system (A: TM — > 
F, b) — then $ is uniquely determined by equation ( \5. jp . ■ 

As for the infinitesimal symmetries, we have a similar situation: equation ( f4.2| ) de- 
termines W = ^4*[y T ], and then equation ( |4.1| ) says TA o V = o A, which means 
[V, A] = 0. So we have: 

Proposition 9 Suppose that the system (A: TM — > F, b) is regular, and let V be a vector 
field in M . Then the following statements are equivalent: 

1. V is an infinitesimal symmetry of the associated implicit system D. 

2. V leaves the dynamical vector field X invariant (\V, X] = 0). 

3. There exists an infinitesimal symmetry (V, W) of the linearly singular system {A: TM 
F, b) — then W is uniquely determined as W = A* [V^] . ■ 
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An important case of a regular system is that of a hamiltonian system (M, u, H), where 
w is a symplectic form on a manifold M and H is a hamiltonian function. This defines a 
linearly singular system (a): TM — ► T*M, d-ff), whose dynamics is ruled by the hamiltonian 
vector field Xh- A diffeomorphism <p preserving Xh is not in general a symmetry of the 
hamiltonian system, since it may not be a canonical transformation (symplectomorphism). 
However, we have shown that it defines a symmetry of the system if considered as a linearly 
singular system. 

5.2 Consistent systems 

Definition A linearly singular system (^4: TM — > F, b) is consistent if A is a surjective 
vector bundle morphism. 

The solutions of the corresponding differential equation are the integral curves of the 
primary vector fields, that is, the sections of the affine bundle D — > M; they can be 
expressed as A D + r, where X is a particular primary vector field and T belongs to Ker A. 
So, the invariance of D can be stated in terms of X Q and Ker A: 

Proposition 10 Let (A: TM — > F, b) be a consistent linearly singular system, and let D 
be the associated implicit system. 

1. A diffeomorphism ip: M — ► M is a symmetry of D iff Ker A is invariant by (p and, 
for any primary vector field X Q , ip*[X ] — X D is in Ker A. 

2. A vector field V in M is an infinitesimal symmetry of D iff, for any vector field T in 
Ker A, [V,T] is in Ker ^4, and, for any primary vector field X Q , [V,X ] is in Ker A. 

Proof. The first statement follows from the equivalence between assertions 1 and 2 in 
theorem [l|. The second statement is the infinitesimal version of the first one. ■ 

Note that a consistent system is locally solvable, so the preceding conditions charac- 
terise the dynamic symmetries of the linearly singular differential equation. 

Let us remark that the definition of a consistent system could be slightly more general. 
If the consistency condition for the linear equation A x -u x = b(x) holds at every x G M 
(i.e., M\ = M), then the image of b is contained in Im^4 C F. If this is a subbundle, we 
could safely substitute Iulj4 for F in the linearly singular system, thus obtaining what we 
have called a consistent system. A similar remark can be applied to regular system. 

5.3 Symmetries at the end of the constraint algorithm 

Consider a linearly singular system (A:TM — > F,b). If the constraint algorithm as ex- 
plained in section 2 can be performed on it (in the sense that the regularity assumption is 
satisfied at each step of the algorithm), then the final dynamics is that of a consistent lin- 
early singular system, (Af. TMf — > ImAf, bf), so the preceding proposition may be directly 
applied to it. 
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6 An example: the associated presymplectic system 

In some problems of control theory, equations of the type A(x)x = b(x,u), where u 
represents the control, play a relevant role. Ibort noted — see [ pi 0C| ] — that, from a linearly 
singular system (A:TM — ► F, b), one can define a presymplectic system (F*,u),H) on the 
total space of the dual vector bundle it*: F* — > M. This is as follows. 

If 9m is the canonical 1-form and ujm = —d&M is the canonical symplectic form of 
T*M, one can use the transpose map *A:F* -> T*M of A to define forms on F* by 
pull-back: 

9 = t A*[e M ], LO = t A*[u M ]. (6.1) 
In a similar way one can use the section b: M — > i 7 to define a linear function ff : i 7 * — > R: 

fT(a s ) = (a a ,,6(a;)>. (6.2) 

One can study the relations between both systems. For instance, each solution of the 
linearly singular equation is in correspondence with a family of solutions of the equation of 
motion of the presymplectic system. We shall limit ourselves to study the relation between 
the symmetries of both systems. 

Consider a vector bundle automorphism (99, <3?) of it: F — > M. The contragradient map 
<J> V = i s a vector bundle automorphism of F*, with base map if. 

In analogy to (|2.8j) , <J> V transforms the 1-form and the hamiltonian function: 

$ v [0] = (T$ v ) v o o (* v ) _1 , = # (* V ) -1 - 

Note that H = (Id, b o ir*). A computation shows that 

^[H] = (ld,^[b]o7T*) (6.3) 

which proves that if b is ^-invariant then H is $ v -invariant; the converse is also true, since 
b is determined by H. 

In a similar way, consider for the sake of simplicity the 1-form 9 as a linear function 
9: TF* — > R — we use the same letter to not overload the notations. Then note that 
9 = {rp*,A o Ttt*). Another computation shows that 

$y[e] = {T F .,$y[A]oTn*). (6.4) 

As before, this proves that A is ^-invariant iff 9 is <3> v -invariant. 

As a conclusion, we have: 
Proposition 11 Consider a linearly singular system (A: TM — ► F,b), and the associated 
presymplectic system (F* , — d9, H). Let (</?, be a vector bundle automorphism of F — > 
M. Then: 

1. b is ^-invariant iff H is Q y -invariant. 

2. A is ^-invariant iff 9 is Q y -invariant. 



So if (</?, <3>) is a symmetry of the linearly singular system then <I> V is a symmetry of the 
presymplectic system. ■ 
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7 Conclusions 

In this paper we have studied the symmetries of linearly singular differential equations. 
To do this, we have intended to clarify what a "symmetry" is for a differential equation, 
and we have found that there are several legitimate approaches to this concept. Some 
approaches rely on the geometry of the data defining the differential equation ( "geometric 
symmetries" ) , whereas others simply characterise the symmetry as a transformation that 
maps solutions into solutions ("dynamic symmetries"). As for the type of symmetries, we 
have considered only point symmetries, which are those defined by diffeomorphisms of the 
configuration space (or vector fields, in the infinitesimal case). 

After a general introduction to linearly singular differential equations, we deal with 
two different concepts of geometric symmetry. The first one, the symmetry of an implicit 
system, has a general applicability and relies on the invariance of the subset D C TM that 
defines, implicitly, the differential equation. The second one, the symmetry of a linearly 
singular system, is specific of the systems discussed in this paper. For them, we show that 
both concepts of symmetry are essentially equivalent. We prove it for a general, finite, 
transformation, and also for the case of flows generated by infinitesimal transformations. 

Under appropriate regularity assumptions, the concept of dynamic symmetry of the 
implicit differential equation defined by D C TM is equivalent to the concept of geometric 
symmetry of an appropriate locally solvable system D' C TM'; this also holds for the 
linearly singular case, by means of the constraint algorithm. We show that, under appro- 
priate regularity assumptions, the most general dynamic symmetry for a linearly singular 
equation (on a manifold M) can be locally realised as a geometric symmetry of a linearly 
singular system (on the final constraint manifold Mf). 

As for the tools needed to deal with the infinitesimal symmetries, in the appendix 
we have performed a careful study of infinitesimal transformations and the invariance of 
geometric structures under the action of flows. This has been applied to infinitesimal 
automorphisms of vector bundles, and may be useful to deal with other problems about 
infinitesimal invariance. 



Appendix: calculus with infinitesimal variations 

General aspects 

To deal with infinitesimal symmetries it will be convenient to perform a more general 
study. Let us consider a map /:R x M — > N. We will use the notation f e {x) = f(e,x), 
and we can interpret the maps f e :M^Nas& "variation" of the map f = /(0, — ). The 
corresponding "infinitesimal variation" of / is the map wj: M — > TN defined as 

d 

/; (A.i) 

£ = 



w/(x) = ¥£ 



in other words, it is w/(x) = f'(0,x), where 

/' = T/oJ?:RxM-> TN 
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is the e-derivative of / — here E is the unit vector field on R interpreted as to be on 
R x M. 

Note that wj is a vector field along f a . (Interpreting / as a path in the infinite- 
dimensional manifold of maps from M to N, wj is its tangent vector at e = 0.) It is clear 
that if f £ does not depend on e, then wj = (the converse is obviously not true). 

Given a number c, if we define g(e,x) = f(ce,x) it is easily seen that w g = cw/, and 
in particular the infinitesimal variation of f(—e,x) is — w/. 

Now let us consider another variation g: R x N — > P, and construct the compositions 
g £ o f £ ; they define a variation g • /: R x M — > P: 

(^./)(e,ar) = 5 (e,/(e,x)). (A.2) 

A direct computation in coordinates shows that 

w 9»/ = w 9 fo + T 9o oWf. (A.3) 
Let us show these objects in a diagram: 




M N P 

This result is immediately extended to the composition of three (or more) variations. 

Another immediate consequence is the following one. Suppose that the f e :M^N 
are diffeomorphisms, and set g £ = ff 1 . Then we have 

w g = -Tg o\v f og . 

Finally, in a similar way one can prove that, if h £ = (f £ , f' £ ) then = (wj, wj/), with 
the usual identification T(N x N') = TN x TN'. 



Transformation of maps 

Consider variations f £ : M — ► M' and g £ : N — > A', and a map /i': M' N'. If the g> £ are 
diffeomorphisms, we can construct a family of maps h £ : M — > as 

^ = 57 1 ° ^ ° Ai 

for which one can easily compute Wh- 
in the important particular case where M = M', N = N', and / and g are variations 
of the identity, obviously one has that h = b! and 

w/, = -w s o h + Th Q o w/. (A.4) 

Now let us take fibre bundles tt:M^B and p:N^C. Suppose that the couple 
g £ :M^N and f £ :B — > (7 is a variation of a bundle morphism (g ,fo)- If (9e,fs) are 
bundle morphisms (i.e., p o g £ = f £ o ir) then 

TpoW 9 = Wjo7T. (A. 5) 
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One can consider more particularly the case where g and / are variations of the identity, 
and ask for instance whether a section a: B — > M is transformed into sections, or is left 
invariant, etc. 



Variations defined in terms of flows 



Note that, with slight complications, the preceding definitions and results could be applied 
to variations / defined only on an open subset D C R X M containing {0} x M. This 
remark is especially relevant for what follows. 

Let us consider a vector field X on M, and denote its flow by Fx, so that the integral 



curve with initial condition x is e 



Let us assume for simplicity that X is 



complete, that is, the domain of Fx is R x M; otherwise we would apply the preceding 
remark. It is clear from the definition of the flow that 



F' x =XoF x , 



(A.6) 



F9 



X. 



■ x = Id, w Fx 

Note that applying ( [A. 6 ) one can compute the e-derivative of the expression F £ (F~ £ (x)) = 
x and obtain 

T x F x .X(x)=X(F x (x)), (A.7) 

which indeed tells that X is invariant under its flow. 

Now let us consider a map h Q : M — > N, and vector fields X on M and Y on N. We 
can use their flows to transform h Q as 

h £ (x) =h(e,x) =Fy £ (h (F x (x))). 

Applying ( |A.4| ) to this composition we obtain 

= Th D o X — Y o h D . 

Note therefore that w/j is zero iff the vector fields X and Y are /i D -related. 



(Ai 



(AS 



In [KMS9S] such infinitesimal variations are studied under the name of "generalised 
Lie derivatives", and for instance w/j in ( |A.9| ) is denoted by C^ X Y)h . 
Let us show all these objects in a diagram: 



M 



N 



TM 



TN 



Fx 



M 



h 



N 




Proposition 12 With the previous notations, h Q is invariant under the action of the 
couple of flows iffwh = 0. 

Proof. The direct implication is a trivial consequence of the definition of Wft. For the 
converse, to show that h(e,x) = h(0,x) we will prove that the e-derivative of h, h'(e,x), 
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vanishes at any e (not only at e = 0). This will be a consequence of the following equation: 

h'(e,x) = T ho{F e x(x)) Fy e -(Th o X -Y o h )(F e x (x)). (A.10) 

To prove it, let us introduce some notation. First, we write F and G the flows of X and Y. 

The tangent map of F applied to a couple of vectors ((e,r),u x .) € T £ R x T X M = 
Tu x \(R x M) can be written 

TF((e,r),u x ) = F\e,x)t + T x F £ -u x . (A.ll) 

Changing some terms in ( |A.7| ), we also have 

Y(y) = T GHy) G- £ -Y(G e (y)). (A.12) 

Now let us proceed to compute h! = Th o E. Let us write h = G o h o F, where 
F(s,x) = (e, F(e, x)), h Q (e,x) = (e,h (x)), and G(e,x) = G(—e,x). The chain rule 
applied to these maps, as well as ( |A.11| ) and ( |A.6| ), yield 

ti(e,x) = -Y(G- £ (h (F £ (x))))+T ho{Fs(x)) G- e -T F . (x) h -X(F F -(x)). 
Application of ( |A.12[ ) converts this equation into 

h'(e,x) = -T fto(F£(;E)) G- £ -y(/ lo (F £ (x))) + T fco(F£(;E)) G- £ -T F£(;E) /i -A(F £ ( a ;)). 
= T ho{F e ix)) G- £ -(T Fe{x) h -X(F £ (x)) - (Yoh )(F £ (x))) 
= TG~ £ °(Th oX -Y oh )(F £ (x)), 
= (TG- £ ow h )(F E (x)), 

which is ( |A.1C| ). ■ 



Infinitesimal vector bundle automorphisms 

Let 7r: E — > M be a fibre bundle. An infinitesimal automorphism is a couple of vector 
fields (X,Y) on M and E such that their flows (F £ ,G £ ) are fibre bundle isomorphisms. 

Proposition 13 With the preceding notations, (X, Y) is an infinitesimal automorphism 
iffY is projectable to X. 

Proof. The condition of being morphism is ttoG £ = F £ on, that is to say, ir = F~ £ oiroG £ . 
According to proposition the right-hand side is e-invariant iff Ttt o Y = X o tt, which is 
the condition of projectability. ■ 

Note that under these conditions the vector field X is determined by Y, so we can as 
well say that Y is an infinitesimal automorphism of the fibre bundle. 

From now on let us suppose that tt: E — » M is a vector bundle. A couple of vector fields 
(X,Y) is an infinitesimal vector bundle automorphism if their flows (F £ ,G £ ) are vector 
bundle isomorphisms. Of course Y is projectable to X. Y is said to be a linear vector 
field if (X, Y) is a morphism between the vector bundles tt: E — > M and T-zr: TE — > TM: 

E TE 



7T 



Tvr 



M — TM 
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Proposition 14 Let Y be a ir-projectable vector field on E. Then Y is an infinitesimal 
vector bundle automorphism iff Y is a linear vector field. 

Proof. First remember that a smooth function defined on a vector space is homogeneous 
of degree one iff it is linear. 

So, denoting by my.E — > E the multiplication by A, the fact that G £ is linear is 
equivalent to G e om\ = m\oG e for each A. This can be expressed also as m\ = G~ £ om\oG £ , 
and applying proposition |l2| again, this holds for each e iff Tm^ oY = Y o m\. Written in 
other terms, we have X-Y(u) = Y(Au), where the first product is meant to be with respect 
to the vector bundle structure Ttt: TE — > TM. So on each fibre Y X :E X — > T X ( X )E is a 
homogeneous smooth map, and therefore linear. ■ 



A different proof of this result, using coordinates, can be found in | KMS 93 |. It is also 
interesting to recall that a projectable vector field Y is linear iff it is invariant under the 
action of the Liouville's vector field, [As, Y] = 0. 
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